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Abstract 

In this paper wc present a new mathematical model for the density regions where a specific spectral 
line and its SACs/DACs are created in the Oe and Be stellar atmospheres. In the calculations of final 
spectral line function we consider that the main reasons of the line broadening are the rotation of the 
density regions creating the spectral line and its DACs/SACs, as well as the random motions of the ions. 
This line function is able to reproduce the spectral feature and it enables us to calculate some important 
physical parameters, such as the rotational, the radial and the random velocities, the Full Width at Half 
Maximum, the Gaussian deviation, the optical depth, the column density and the absorbed or emitted 
energy. Additionally, we can calculate the percentage of the contribution of the rotational velocity and the 
ions' random motions of the DACs/SACs regions to the line broadening. Finally, we present two tests and 
three short applications of the proposed model. 
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1. Introduction 

When we study the UV lines of hot emission stars we 
have to deal with two problems: (i) The presence of a 
very complex structure of many spectral lines in the UV 
region, such as the resonance lines of Si IV, C IV, N 
V, Mg II and the N IV spectral line, (ii) The pres- 
ence of Discrete Absorption Components (DACs, Bates 



& Halliwell (1986)). 

DACs are discrete but not unknown absorption spectral 
lines. They are spectral lines of the same ion and the same 
wavelength as a main spectral line, shifted at different 
AA, as they are created in different density regions which 
rotate and move radially with different velocities (Danezis 
1983; Danezis 1987; Danezis et al. 2003; Lyratzi et al. 
2007). DACs are lines, easily observed, in the case that 
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the regions that give rise to such hnes, rotate with low 
velocities and move radially with high velocities. 

Many suggestions have been made in order to explain 
the DACs phenomenon. Most of researchers have sug- 
gested mechanisms that allow the existence of structures 
which cover all or a significant part of the stellar disk, such 
as shells, blobs or puffs (Underhill 1975; Henrichs 1984; 
Underhill & Fahey 1984; Bates & Halliwell 1986; Grady 
et al. 1987; Lamers et al. 1988; Waldron, Klein, & Altner 
1992; Waldron, Klein, & Altner 1994; Cranmer & Owocki 
1996; Rivinius et al. 1997; Kaper et al. 1996; Kaper 
et al. 1997; Kaper et al. 1999; Markova 2000), inter- 
action of fast and slow wind components, Corotation 
Interaction Regions (CIRs), structures due to magnetic 
fields or spiral streams as a result of the stellar rota- 
tion (Underhill & Fahey 1984; MuUan 1984a; MuUan 
1984b; MuUan 1986; Prinja & Howarth 1988; Cranmer 
& Owocki 1996; FuUerton et al. 1997; Kaper et al. 
1996; Kaper et al. 1997; Kaper et al. 1999; Cranmer, 
Smith, & Robinson 2000). Though we do not know yet 
the mechanism responsible for the formation of such struc- 
tures, it is positive that DACs result from independent 
high density regions in the stellar environment. 

An important question is whether there is a connec- 
tion between the observed complex structure of the above 
spectral lines and the presence of DACs. A possible an- 
swer is that if the regions that create the DACs rotate 
with large velocities and move radially with small ve- 
locities, the produced lines have large widths and small 
shifts. As a result, they are blended among themselves 
as well as with the main spectral line and thus they 
are not discrete. In such a case the name Discrete 
Absorption Components is inappropriate and we use 
only the name Satellite Absorption Components (SACs) 
(Lyratzi & Danezis 2004; Danezis et al. 2005; Danezis et 
al. 2006; Nikolaidis et al. 2006a; Nikolaidis et al. 2006b). 
A very important question is whether it is possible that 
the existence of SACs may be responsible for the complex 
structure of the observed spectral feature. 

As it is clear, for a future study of the mechanisms 
responsible for the creation of the density regions which 
produce the complex profiles of the above spectral lines, as 
well as for the study of their structure and evolution, we 
need to calculate the values of some physical parameters 
of these regions. 

In order to calculate these parameters, it is necessary 
to construct a line function, based on the idea of DACs 
and SACs phenomena able to reproduce, in the best way, 
the observed spectral feature. 

Danezis et al. (2003) presented such a line function in 
the case that the reason of the line broadening is only the 
rotation of the regions which create the observed spectral 
lines and their components and when these regions present 
spherical symmetry around their own center or the center 
of the star (see also Lyratzi et al. (2007)). With this line 
function we calculate some important physical parame- 
ters, such as the rotational and the radial velocities, the 
optical depth, the column density (Danezis et al. 2005) 
and the absorbed/emitted energy. But, in some cases the 



chaotic motion of emitters/absorbers inside the dense lay- 
ers can be significant (see e.g. Danezis et al. (2006)) and 
the estimated rotation might be overestimated. Therefore, 
the modification of the previously given model (Danezis et 
al. 2003) is needed, in sense that the possible contribution 
of random velocities to the line broadening can be taken 
into account. 

Here we present a new model which includes also the 
random motions of the ions, since in the case of hot emis- 
sion stars one can expect a significant contribution of ran- 
dom motion of the emitters/absorbers to the line profile. 
Based on this idea, besides the above mentioned physical 
parameters of the regions that create the complex spectral 
lines, we can calculate the mean random velocity of the 
ions, but also the percentage of the contribution of the 
rotational velocity and emitter/absorber random motions 
of the DACs/SACs regions to the line broadening, as well 
as the Gaussian deviation for each DAC/SAC profile. 

The paper is organized as follows: In §2 we give the 
description of the model, in §3 we discuss the model, in 
§4 we give two tests of the model, in §5 we present ap- 
plications of the model and finally, in §6 we outline our 
conclusions. 

2. Description of the model - The Line Function 

As was already mentioned above, Danezis et al. (2003) 
presented a line function which is able to reproduce accu- 
rately the observed spectral lines and their SACs/DACs 
in the same time. The proposed line function is: 
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where: 

Fq (A): the initial radiation flow, 

Li, Lej, Lg-. are the distribution functions of the absorp- 
tion coefficients k\i, kx^j, k\g respectively. Each L de- 
pends on the values of the apparent rotational velocity as 
well as of the radial velocity of the density shell, which 
forms the spectral line (K-ot, Vrad) 
^f. is the optical depth in the center of the line, 
Sxej : the source function, which, at the moment when the 
spectrum is taken, is constant. 

In Eq. 1, the functions Sxej {'^ - e^^"^^'^) , 

Q-Lg£,g g^]-g ^YiB distribution fimctions of each satellite com- 
ponent and we can replace them with a known distribu- 
tion function (Gauss, Lorcntz, Voigt or disk model). An 
important fact is that in the calculation of F{\) wc can 
include different geometries (in the calculation of L) of 
the absorbing or emitting independent density layers of 
matter. 

The decision on the geometry is essential for the cal- 
culation of the distribution function that we use for each 
component, i.e. for different geometries we have different 
line shapes, representing the considered SACs. 

Eq. 1 gives the function of the complex profile of a spec- 
tral line, which presents SACs or DACs. This means that 
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Eq. 1 is able to reproduce not only the main spectral line, 
but its SACs as well. 

The main hypotheses when we constructed the Rotation 
distribution function were that the line width AA is only 
due to the rotation of the regions which create the ob- 
served spectral lines and their components and that these 
regions present spherical symmetry around their own cen- 
ter or the center of the star (see also Lyratzi et al. 
(2007)). Consequently, the random velocities of the emit- 
ters/absorbers in the density region are assumed to be 
negligible, i.e. they do not significantly contribute to the 
line profile. Here we consider that the random velocities 
may have a significant contribution to the distribution 
function L (see Appendix 1) In this case the final form 
of the distribution function L is given as: 




where: 

Ao is the transition wavelength of a spectral line that arises 
from a specific point of the equator of a spherical density 
region that produces one satellite component, z = 
{Vrot is the rotational velocity of the specific point) and 

erf (x) — - J e~" du is the function that describes the 



Gaussian error distribution. 

3. Discussion and way of application of the pro- 
posed model 

Introducing the previous final function of a complex 
spectral line (Eq. 1), in combination with the distribu- 
tion function L given in Eq. 2, we would like to note the 
following: 

1. The proposed line function (Eq. 1) can be used for 
any number of absorbing or emitting regions. This means 
that it can also be used in the simple case that i = 1 and 
J = or z = and j = 1 , meaning when we deal with sim- 
ple, classical absorption or emission spectral lines, respec- 
tively. This allows us to calculate all the important phys- 
ical parameters, such as the rotational, the radial and the 
random velocities, the Full Width at Half Maximum, the 
Gaussian deviation, the optical depth, the column density 
(Danezis et al. 2005) and the absorbed or emitted energy, 
for all simple and classical spectral lines in all spectral 
ranges. 

2. For each group of the parameters Vroti, Vradi, Vrandi 
and .fi, the function I\i = e^^^^' reproduces the spectral 
line profile formed by the i density region of matter, mean- 
ing that for each group we have a totally different profile. 
This results to the existence of only one group of Vroti, 
Vradi, Vrandi and S.i giving the best fit of the i component. 
In order to accept as the best fit of the observed spectral 
line, given by the groups (Vrou, Vrad,, Vrand,, 6) of all 
calculated SACs, one has to adhere to all physical criteria 
and techniques, such as: 



i) To make a complete identification of spectral lines in 
the region around the studied spectral line and to have the 
superposition of the spectral region that we study with the 
same region of a classical star of the same spectral type 
and luminosity class, in order to identify the existence of 
spectral lines that blend with the studied ones, as well as 
the existence of SACs. 

ii) The resonance lines as well as all lines originating 
in a particular region should have the same number of 
SACs, depending on the structure of this region, with- 
out influence of ionization stage or ionization potential 
of emitters/absorbers. As a consequence, the respective 
SACs should have similar values of radial and rotational 
velocities. 

iii) The ratio of the optical depths of two resonance lines 
must be the same as the ratio of the respective relative 
intensities. 

3. In order to conclude to the group of the parameters 
which give the best fit, we use the model in the following 
two steps: 

(i) At the first step we consider that the main reason 
of line broadening of the main line and satellite compo- 
nents is the rotation of the region which creates the com- 
ponents of the observed feature and a secondary reason 
is the thermal Doppler broadening. This means that we 
start fitting the line using the maximum Vrot- Then we in- 
clude Doppler broadening, in order to accomplish the best 
fit (Rotation case), (n) At the second step, we consider 
the opposite. In this case the main reason of line broaden- 
ing of the main line and satellite components is supposed 
to be the Doppler broadening and the secondary reason 
is rotation of the region which creates components of the 
observed feature. This means that we start fitting the line 
using the maximal Doppler broadening. Then we include 
Vrot, in order to accomplish the best fit (Doppler case). 

In both cases (Rotation case and Doppler case) we check 
the correct number of satellite components that construct 
the whole line profile. At first we fit using the number of 
components that give the best difference graph between 
the fit and the real spectral line. Then we fit using one 
component less than in the previous fit. The F-test be- 
tween them allows us to take the correct number of satel- 
lite components that construct in the best way the whole 
line profile. The F-Test between these two cases indicates 
the best way to fit spectral lines. When the F-Test can- 
not give definite conclusion on which case we should use, 
we still could obtain information about limits of Vrot and 
a. If the F-Test gives similar values, then the rotation 
case defines the maximal Vrot and the minimal a and the 
Doppler case defines the minimal Vrot and the maximal a. 

4. Profiles of each main spectral line and its SACs 
are fitted by the function e~^'^'in the case of an absorp- 
tion component or Sxej (l — e~^=J^=^) in the case of an 
emission component. These functions produce symmetri- 
cal line profiles. However, most of the spectral lines are 
asymmetric. This fact is interpreted as a systematical 
variation of the apparent radial velocities of the density 
regions where the main spectral line and its SACs are cre- 
ated. In order to approximate those asymmetric profiles 
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we have chosen a classical method. This is the separation 
of the region, which produces the asymmetric profiles of 
the spectral line, into a small number of sub-regions and 
each of them is treated as an independent absorbing shell. 
In this way we can study the variation of the density, the 
radial shift and the apparent rotation as a function of the 
depth in every region that produces a spectral line with an 
asymmetric profile. All mentioned above have to be taken 
into account during the evaluation of our results and one 
should not consider that the evaluated parameters of those 
sub-regions correspond to independent regions of matter, 
which form the main spectral line or its SACs. 

5. We suggest that the width of the blue wing is the re- 
sult of the composition of profiles of the main spectral line 
and its SACs. Thus, the blue wing of each SAC gives the 
apparent rotational velocity of the density shell, in which 
it forms. In order to have measurements with physical 
meaning, we should not calculate the width of the blue 
wing of the observed spectral feature but the width of the 
blue wing of each SAC. 

6. We would like to point out that the final criterion to 
accept or reject a best fit, is the ability of the calculated 
values of the physical parameters to give us a physical 
description of the events developing in the regions where 
the spectral lines presenting SACs are created. 

7. In the proposed distribution function an important 
factor is m = . This factor indicates the kind of the 
distribution function that is able to fit in the best way 
each component's profile. 

i) If m ~ 3 we have equivalent contribution of the rota- 
tional and random motions to the line widths. 

ii) If m ~ 500 the line broadening is only an effect of the 
rotational velocity and the random velocity is negligible. 
In this case the profile of the line is the same with the 
theoretical profile deriving from the Rotation distribution 
function. 

iii) Finally if m < 1 the line broadening is only an ef- 
fect of random velocities and the line distribution is the 
Gaussian. 

4. Testing the model 

In order to check the validity of our model we perform 
two tests: 

I) In order to check the above spectral line function, 
we calculated the rotational velocity of He I A 4387.928 
A absorption line for five Be stars, using two methods, 
the classical Fourier analysis and our model. In Fig. 1 we 
present the five He I A 4387.928 A fittings for the stud- 
ied Be stars and the measured rotational velocities with 
both methods. The obtained rotational velocities from 
our model are in good agreement with ones obtained with 
Fourier analysis. 

The values of the rotational velocities, calculated with 
Fourier analysis, some times, may present small differ- 
ences compared to the values calculated with our method, 
as in Fourier analysis the whole broadening of the spec- 
tral lines is assumed to represent the rotational velocity. 



In contrary, our model accepts that a part of this broad- 
ening arises from the random motion of the ions. 

We point out that with our model, apart from the rota- 
tional velocities, we can also calculate some other param- 
eters, such as the standard Gaussian deviation ct, the ve- 
locity of random motions of emitters/absorbers, radial ve- 
locities of the regions producing studied spectral lines, the 
full width at half maximum (FWHM), the optical depth, 
the column density and the absorbed or emitted energy. 

II) Additional test of our model is to calculate the ran- 
dom velocities of layers that produce the C IV satellite 
components of 20 Oe stars with different rotational veloc- 
ities. 

We analyzed the C IV line profiles of 20 Oe stars 
which spectra were observed with the lUE - satellite 
(lUE Archive Search database^). We examine the com- 
plex structure of the C IV resonance lines (AA 1548.155, 
1550.774 A). Our sample includes the subtypes 04 (one 
star), 06 (four stars), 07 (five stars), 08 (three stars) 
and 09 (seven stars). The values of the photospheric ro- 
tational velocities are taken from the catalogue of Wilson 
(1963) (see also Antoniou et al. (2006)). 

In the composite C IV line profiles we detect two com- 
ponents in 9 stars, three in 7 stars, four in 3 stars and five 
in one star. 

In Fig. 2 we present the C IV resonance lines best fit 
for the star HD 209975. 

In Fig. 3 we present the random velocities (Vrand) of 
each SAC as a function of the photospheric rotational ve- 
locity {Vphot) for all the studied stars. 

As one can see in Fig. 3 the obtained values for ran- 
dom velocities are in accordance with the classical theory; 
the values of the random velocities do not depend on the 
inclination angle of the rotational axis. As the ioniza- 
tion potential of the regions where satellite components 
are created is the same for all studied stars, one can ex- 
pect similar average values of the random velocities for 
each absorbing region (here denoted as a-d component, 
see Fig. 3) for all studied stars, as we obtained by using 
the model. 

Both of the tests, described above, support our new 
approach in investigation of the regions around hot stars, 
implying that besides rotation of the regions, one should 
consider also the random motion of emitters/absorbers. 

5. Applications of the model: Rotation - density 
regions vs. photosphere 

Here, we give some examples of application of the 
model. We apply the model to study the complex struc- 
ture of C IV and Si IV spectral fines in Oe and Be stars, 
as well as the radial and rotational velocities of different 
components, in order to find more precise rotational com- 
ponent. 

^ http;/ /archive. stsci.edu/iuc/search.php 



No. ] A new model for the structure of the DACs and SACs regions in the Oe and Be stellar atmospheres 5 



HD 56139 
He IX 4387.928 A 




HD 57219 
He 1 1 4387.928 A 



Fourier: Vrot = 88 km/s 
GR model: Vrot = 84 km/s 
o = 0.52 




Fourier: Vrot = 90 km/s 
GR model: Vrot = 90 km/s 
o = 1.22 




Fig. 2. The C IV resonance lines (AA 1548.155, 1550.774 A) 
best fit with GR model for the star HD 209975. The compo- 
nents obtained from the best fit are shown bottom. 



HD 75311 
He IX 4387.928 A 




Fourier: Vrot = 257 km/s 
GR model: Vrot = 260 km/s 
<5= 1.15 



HD 77320 




He I X 4387.928 A 






A. 










Fourier: Vrot = 338 km/s 




GR model: Vrot = 338 km/s 




o = 1.3 



HD 201733 
He I X 4387.928 A 



Fourier: Vrot = 351 km/s 
GR model: Vrot = 345 km/s 
<5= 1.5 



Fig. 1. The five Hel A 4387.928 A fittings for the studied Be 
stars with RG model and the measured rotational velocities 
with Fourier analysis and RG model. The differences between 
the observed and the reproduced spectral lines are hard to see, 
as we have accomplished the best fit. 



5.1. C IV density regions of 20 Oe stars 

In this application we use the previously mentioned C 
IV lUE spectra of the second test. We study the relation 
between the ratio Vrot/yphot of the first, second, third 
and fourth detected component with the photospheric ro- 
tational velocity (Vphot)- This ratio indicates how much 
the rotational velocity of the specific C IV layer is higher 
than the apparent rotational velocity of a star (see also 
Antoniou et al. (2006)). In Fig. 4 we present our results. 

In each region and for each component we can conclude 
that there exists an exponential relation between the ratio 
Vrot /Vphot and the photospheric rotational velocity Vphot- 
The maximum ratio Vrot /Vphot varies from 40, for the first 
to 5 for the fourth component (Fig. 2). A possible expla- 
nation of this situation is the inclination of the stellar axis. 

5. 2. Si IV density regions of 27 Be stars 

This study is based on the analysis of 27 Be stellar spec- 
tra taken with the lUE - satellite. We examine the com- 
plex structure of the Si IV resonance lines (AA 1393.755, 
1402.77 A). Our sample includes all the subtypes from BO 
to B8. The values of the photospheric rotational velocities 
are taken from the catalogue of Chauville et al. (2001). 

We found that the Si IV spectral lines consist of three 
components in 7 stars, four in 15 stars and five in 5 stars. 
We study the relation between the ratio Vrot /Vphot of the 
first, second, third fourth and fifth detected component 
with the photospheric rotational velocity {Vphot)- This 
ratio indicates how much the rotational velocity of the 
specific Si IV layer is higher than the apparent rotational 
velocity of the star (see also Lyratzi et al. (2006)). In 
Fig. 5 we present our results. 

The Si IV resonance lines are composed of three four or 
five components. The difference with the case of the C IV 
resonance lines in the spectra of 20 Oe stars is that they 
are composed of two, three or four components. However, 
in both cases, in each region and for each component there 
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Fig. 3. Random velocities (Vrand) of the four SACs as a 
function of the photospheric rotational velocity (V^^ot) for 
all the studied Oo stars. 

exists an exponential relation between the ratio Vrot /Vphot 
and the photospheric rotational velocity Vphot- For the 
satellite components of the Si IV resonance lines, the max- 
imum ratio Vrot /Vphot varies from 19, for the first to 1 for 
the fifth component (Fig. 5). The same phenomenon ap- 
pears in the case of the C IV resonance lines in 20 Oe stars, 
where the maximum ratio Vrot /Vphot varies from 40, for 
the first to 5 for the fourth component. 

5.3. Rotation in photosphere vs. rotation of density re- 
gions 

As one can see in Figs. 4 and 5, there is a good corre- 
lation between the rotational velocities of density regions 
and photosphere. In both cases (Be and Oe stars) for 
rotational velocity of the photosphere V^hot > 200 km/s, 
the rotation of the photosphere correlates with the rota- 
tion of the density regions (Vi-ot ~ const, x V^hot)- On the 



Fig. 4. Ratio VrotVpfiat of the four SACs as a function of the 
photospheric rotational velocity (Vphot) for all the studied Oe 
stars. 

other hand, for Vphot < 200 km/s, the ratio Kot/V'phot de- 
creases with the rotation in the photosphere. This can 
be explained by the inclination effect and the fact that 
the density regions are extensive around the star, i.e. for 
a high inclination angle the projected photospheric rota- 
tional velocity is small, but in a density region (that lies 
extensively around the star), one can detect faster rota- 
tion. 

6. Conclusions 

Here we present a new model for fitting the complex UV 
lines of Be and Oe stars, where we take into account the 
possibility that random motion of ions can significantly 
contribute to the line widths. The proposed model can 
be applied to the spectral photospherical lines as well 
as the UV lines originating in the post-coronal regions. 



No. 



A new model for the structure of the DACs and SACs regions in the Oe and Be stellar atmospheres 



7 



5 
> 



Vrot/Tphot 

a compoueut 



2D - 
15 - 
ID - 



:00 3DD 
I'phot (kur's) 



ViotA'phot 

b coaipoueut 



100 200 300 400 500 

Vphot (kai s ) 



ViotA^phot 

c component 



* • • * 



100 :oo 300 

\ phot (kni s) 



400 



500 



Viot/Vphot 

d coQipone ut 



1.5 



100 :oo 300 

^ phot (kni's) 



400 



500 



ViotA'phot 

e conipone nt 



1.5 
1 

O.S 




100 200 300 400 

Vphot (km s) 



Concerning our work we can give the following conclu- 
sions: 

(i) The proposed model can reproduce accurately the 
complex UV lines of Be and Oe stars. 

(ii) Using the proposed model one can very well separate 
the contribution of the rotational and radial velocities in 
the density region fitting the complex line profiles. 

(iii) The proposed model gives the opportunity to in- 
vestigate physical parameters of the regions creating the 
UV complex line profiles. 

(iv) The proposed model allows us to use the UV pho- 
tospheric lines, in order to determine the photospheric 
rotation. 

At the end let us point out that in spite of the fact that 
today exists a number of models which are able to repro- 
duce spectra from stellar atmospheres, there is a problem 
to find appropriate model that can fit the complex UV line 
profiles which are created not only in the photosphere, but 
also in the density regions. On the other hand, the pro- 
posed model is able to provide the information about the 
physical parameters of density regions as well as rotation 
of the photosphere. We hope that the proposed model will 
be useful first of all to have first impression about physics 
of density layers. 
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Fig. 5. Ratio Vrot/Vpfiot of the five SACs as a function of the 
photospheric rotational velocity (Vphot) for the studied Be 
stars. 
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Appendix 1. Including the random motion in the calculation of L - The Gauss- Rotation model (GR 
model) 



We consider a spherical shell and a point Ai in its equator. If the laboratory wavelength of a spectral line that arises 
from Aj is A;o6, the observed wavelength is: Aq = Xiab ^ ^Kad- If the spherical density region rotates, we observe 
a displacement AArot and the new wavelength of the center of the line is A^ = Aq ± AArot where AArot = Xozsinip, 
z = iiiit, Vrot is the rotational velocity of the point A^. This means that A^ = Aq ± X^zsaup = Aq (1 ± zsvtn^) and if 
— ^ < (/3 < ^ then A; = Ao (1 — zsaup). If we consider that the spectral line profile is a Gaussian distribution, then 

/(A)- 1 - 



27rcr 



where k is the mean value of the distribution and in the case of the line profile it indicates the 



center of the spectral line that arises from A^. This means that: 



/(A) = 



1 



^ — ^0 (1 — ^ siny) "I 



1 



The distribution function for all the semi-equator is: 



/i(A) 



1 



[A-A0(l-zsiny)l^ 



27rcr 



' cosifdtp 



(Al) 



If siniy9 = X then dx ~ cosipdip, — 1 < a; < 1, and Eq. Al takes the form 



/i(A) 



If we set 



1 



1 a^/2TT 

A — Ao{l— ^a::) 

A-Apd- 



/i(A) 



1 



XozJn 



2cr2 

we have 



A-Ao(l + z) 

CT\/2 



We consider the function erf (x) = §/ e~" du. It is the known function which describes the Gaussian error distri- 



bution. 

If we take into account this function, /i(A) takes the form 



Ii(A) 



Ii(A) 



1 



XqzJtt 



Xoz^/^T 



A-Ao(l~z) 

tT>/2 



A-AqCI + z) 
!tV2 



du - 



du 



TT fX-Xo{l-z)\ TT fX-Xo{l + z) 

—erf ^ erf ^ 



Thus, we finally have 



Ii(A) 



2Xoz 



,,A-Ao(l-2)\ /A-Ao(l + 2) 

erj ■= — erj ' 



(7^2 



aV2 



and the distribution function from the semi-spherical region is: 



I/iriai (A) — 



2Xoz 



erf 



A-Aq ^ 

aV2 (tV2 



Ao2 a\ „/A-Ao Xqz 

cosa —erf | ^ ^cos( 



aV2 G\f2 



cos QdO 



(A2) 



(Method Simpson) 

In Eq. A2, from Aq we calculate the value of the radial velocity (K-ad), from z we calculate the rotational velocity 
(K-ot) and from cr we calculate the random velocity {Vrand)- 

This distribution function \ final (A) has the same form with the distribution function of the absorption coefficient 
L and may replace it in the line functions e~^^ or S\f^ (l — e~^=^^=^), in the case when the line broadening is an effect 
of both the rotational velocity of the density region as well as the random velocities of the ions. This means that now 
we have a new distribution function to fit each satellite component of a complex line profile that presents DACs or 
SACs. We name this function Gauss- Rotation distribution function (GR distribution function). 
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